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Recently a new optimal control modification has been introduced that can achieve robust adaptation with a 
large adaptive gain without incurring high-frequency oscillations as with the standard model-reference adap- 
tive control. This modification is based on an optimal control formulation to minimize the f£j_ norm of the 
tracking error. The optimal control modification adaptive iaw results in a stable adaptation in the presence of 
a large adaptive gain. This study examines the optimal control modification adaptive law in the context of a 
system with a time scale separation resulting from a fast plant with a slow actuator. A singular perturbation 
analysis is performed to derive a modification to the adaptive law by transforming the original system into a 
reduced-order system in slow time. The model matching conditions in the transformed time coordinate results 
in increase in the feedback gain and modification of the adaptive iaw. 


I. Introduction 

In recent years, adaptive control has been receiving a significant amount of attention. There has been a steady 
increase in the number of adaptive control applications in a wide range of settings such as aerospace, robotics, process 
control, etc. 1,2 The ability to accommodate system uncertainties and to improve fault tolerance of a control system 
is a major advantage of adaptive control. Nonetheless, adaptive control still faces significant challenges in providing 
robustness in the presence of unmodeled dynamics and parametric uncertainties. 3 The ability for an adaptive control 
algorithm to modify a pre-existing control design is considered a strength and at the same time a weakness. The 
crash of the X-15 aircraft in 1967 4 serves as a reminder that adaptive control is still viewed with some misgivings 
despite enormous advances ever since. Over the past several years, various model-reference adaptive control (MRAC) 
methods have been investigated. 5-9 

In the conventional MRAC framework, the tracking error is generally inversely proportional to the magnitude of 
the adaptive gain. However, a large adaptive gain can lead to high-frequency oscillations which can excite unmodeled 
dynamics that could adversely affect the stability of an MRAC law. 10 Various modifications were developed to increase 
robustness of MRAC by adding damping to the adaptive law to reduce high-frequency oscillations. Two well-known 
modifications in adaptive control are the cr-modification 1 1 and £i- modification. 12 These modifications have been 
used extensively in adaptive control. Recently, a new modification has been introduced that is based on an optimal 
control formulation to minimize the Aft-norm of the tracking error. 13 The optimality condition results in a damping 
term in the adaptive law proportional to persistent excitation. The optimal control modification has been shown to be 
able to achieve fast adaptation with a large adaptive gain without compromising stability robustness while preserving 
tracking performance. This study extends the development of the optimal control modification adaptive law to the case 
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when there exists a time-scale separation between a fast plant and a slow actuator which prevents the plant to follow 
a reference model even in the presence of adaptive control. A singular perturbation approach is used to separate the 
time scales of the plant and actuators and then modify the optimal control modification adaptive law to account for 
the slow actuator in the singularly perturbed system. The singular perturbation approach transform the original system 
into a reduced-order system in slow time. The model matching condition is applied to the reduced-order system 
and the reference model in the transformed slow time coordinate that results in changes in the actuator command to 
accommodate the slow actuator dynamics. The resulting control signal can then track the reference model better than 
if the actuator command is not modified. 

II. Singularly Perturbed Systems with Slow Actuators 

A direct MRAC problem is posed as follows: 

Given a nonlinear plant as 

x = Ax + B [m + ©* t <!>(x) + v(f)l (1) 

where x(t) : [0,°°) — > R" is a state vector, u ( t ) : [0,°°) — > R" is a control vector, A G R" x " and B € R' ,x " are known 
matrices such that the pair (A . B) is controllable and furthermore A is Hurwitz, ©* G R px " is an unknown constant 
weight matrix, and <t> (x) : R' 1 — > R p is a known bounded bounded basis function and is at least piecewise smooth in x, 
and v (. t ) : [0, °°) — » R" is a small unknown bounded disturbance and is differentiable with || v (?) || < vo G R for all t. 
The controller u (?) is subject to linear dynamics 

u = eA (u — u c ) (2) 

where u c (?) : [0,°°) — > R" is an actuator command vector, e is a positive constant, and A G R' ,x " is a known Hurwitz 
matrix. 

The objective is to design the controller u ( ? ) that enables the plant to follow a reference model 

Am = A m x /n -f- B m r (3) 

where A m G R" x " is Hurwitz and known, B m G R” x " is also known, and r (?) : [0,°°) — > R" G Jzfoo is a command vector 
with r G J§foo. 

If the actuator dynamics are sufficiently fast relative to the reference model dynamics, that is, e ||A|| ||A m ||, then 

the effect of actuator dynamics may be negligible. Then we design a controller u ( t ) to follow an actuator command as 

u c = K x x + K r r - u ad (4) 

where K x G R” x ” and K, G R' ,x " are known nominal gain matrices, and u ad G M" is a direct adaptive signal. 

Defining the tracking error as e = x m — x, then the tracking error equation becomes 

e = x m — x = A m x m + B m r — A m x + A m x — Ax — BK x x — BK r r + B u ad — © 1 <I> (x) — 8 (t) (5) 

Assuming the gain matrices K x and K, can be chosen to satisfy the model matching conditions A + BK X = A„, and 
BK r = B m , and the adaptive signal u ad chosen as 

Uad = (x) (6) 

where © is an estimated weight matrix and © = © — ©* is a a weight variation, then 

e=A m e + B © T <I>(x) — v(f) (7) 

The unknown weight matrix © can be estimated by a standard cr-modification model-reference adaptive control as 

© = -r ( <S>e J PB + (7©) (8) 

In this study, we are interested in the case of slow actuator dynamics when g < 1 is a small parameter and 
e || A|| <C ||A||. Then x(f) is a fast state and u (?) is a slow state. To decouple the fast and slow states, we perform 
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a time-scale separation by applying the singular perturbation method. Toward that end, we consider a slow time 
transformation 

X = £t (9) 

where T is a slow time variable. 

Then, the plant and actuator models are transformed into a singularly perturbed system as 


dx 

e — =Ax + B 
dx 


M + ©* T <t>(x)+v(/) 


( 10 ) 


du 

— = A(u-u c ) 

dx 


( 11 ) 


The Tikhonov’s theorem can be used to approximate the solution of the singularly perturbed system with the 
solution of a “reduced-order” system by setting e = 0. 14 Then, x(u, e) is on a fast manifold. Thus, the reduced-order 
system is given by 

fl" 1 Ar 0 + Mo + ©* T< l > (*o)+v(^ =MoTv^j +/(*o) =0 (12) 


du 0 . , , 

— — = A (mo -M e) 
a t 


(13) 


where xo and «o are the “outer” solution of the singularly perturbed system. 

The term “outer” is in connection with the concept of “boundary layer” or “inner” and “outer” solutions which 
have the origin in boundary layer theory due to Prandtl. The “boundary layer” solution for the singularly perturbed 
system is defined by 


Xi =A(xo+Xi)+B (mq + Ui) + ©* ! <f > (xq + Xi) + v ( ^ 


—Axq — B mo + ©* t< I>(xo) + v^-^ 

= Axj + B L + ©* T <1> (x 0 + x,') - ©* T <I> (x 0 )l (14) 

Ui = £A(mo + Mi — Me) — EAmo = £A(m; — M e ) (15) 

The solution of the original system can be obtained by a matched asymptotic expansion method applied to both 
the inner and outer solutions. 15 The outer solution is in fact the asymptotic solution of the original system as t — > 

The algebraic solution of Eq. (12) can be expressed in general as 


*o=*(«o+v(i)) = -r‘(»o+v(i)) 


(16) 


assuming / 1 exists. 

Differentiating Eq. (16) with respect to the slow time variable and then substituting the actuator model into the 
result yield 

dx o _ dg dup dg d8 _ dg _ dgc/v 
dx duo dx 98 dx duo 0 c dv dx 


From Eq. (12), we have 

mo = — B~^Axq — ©* *£> (x'o ) — 17 ^ ~ ^ 

Hence, we obtain the following reduced-order plant model constrained by the slow actuator dynamics 

— ZT^Axo — ©* T< f>(x 0 ) — v u c 


(17) 

(18) 


dx o _ dg 
dx 9uq 


dg_dv_ 
dv dx 


From Eq. (12), we also have 


fi -' A - l -©* 76 

b-'a + ©* t< 


Thus 


dx o _ dg 
duo duo 


d® (x 0 y 

dx o _ 

dx o 

duo 

■d<P(xo) 


dx o 

dv 

B- l A + ®* Td<t>{Xo) ] 
dx o 


(19) 

( 20 ) 
( 21 ) 

( 22 ) 
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dxo 

dv 


dg 

dv 


's-‘A+e- T ‘ M,(j(l) 


dx o 


Therefore 


Let 


jg _ <5g 

<9 mq <9v 


a, = b-'aab~ 1 a 


(23) 

(24) 

(25) 


and if A s is Hurwitz, then the Tikhonov’s theorem guarantees that the reduced solution with e > 0 converge to the 
solution of the original system with e = 0 as £ — > 0. 14 

Note that Eq. (19) satisfies the outer solution of the nonlinear plant model and the actuator dynamics. Because of 
the slow actuators, the time scale of the response of the plant cannot exceed that of the actuators. Thus, if the reference 
model is faster than the actuator model, the tracking error cannot be guaranteed to be small even with adaptive control 
due to the model mismatch. A possible solution is to revise the reference model to match the actuator-constrained 
plant model, or alternatively to re-design the actuator command to reduce the tracking error. 

In this study, we will consider asymptotic solution of the singularly perturbed system. In effect, the inner solution 
is neglected so that 

x ss xq (26) 


u sa mo 


In slow time, the reference model is expressed as 


/7v 1 

= -(A m x m + B m r ) 

dx e 


(27) 


(28) 


Note that since dg/duo contains the uncertainty, the control design is quite complicated. In order to simplified the 
solution, we make an assumption that the uncertainty term is small. That is 


© 


*T 


d<t> (x) 


dx 


Then, using the matrix inversion lemma, we obtain 
n -1 


b~'a+ 


dx 


= A- l B-A~ 1 B 


© 


*T 


d<i>(x) 

dx 


< llB-UI 


+ {B~ l A\ 


(29) 


-l 


a-'b 


t A~ l B -A- 1 B®* t <1<P ^ A~ l B (30) 
dx 


Then, we make the following choice for the actuator command signal 

u c = K x x + K r r u a d 

where 

K x = A - 1 B~ 4 1 A —A m — B~ l A 


(31) 


(32) 


K r =A~ l B~ l A-B„ 


(33) 


Comparing this controller with the controller when actuator dynamics are fast, the increase in the control gain is 
estimated as 

ll^ll IWI 


ll*; 


e||A|| 


(34) 


where K* is the control gain for fast actuator dynamics. 

The closed-loop singularly perturbed system now becomes 


dx 

dx 


B-'A + Q 


,*T 


d<$> (x) 
dx 


1 1 

B~ l A— (. A m x+B m r ) - 
e 


b-'a+q * 7 


d<i>(x ) 
dx 


u ad -®* J <S>{x) 


dg 

dv 


dv 

dx 


-Av 


(?)' 


(35) 
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Using the result of the matrix inversion lemma, we get 


dx 

dx 


i-a 1 b@* 

dx 


1 


(A m x B m v ) 


b-'a + 0 *t^W 

dx 


-l 


Uad -®* T <t>(x ) 


+ 


dg 

dv 


dv 

dx 


-Av 


(?) 


(36) 


Then, the adaptive signal u ad can be designed to keep the following expression small by a judicious choice of new 
basis function <f>i (x. r ) that spans the unknown parameter space ©j such that 


-A- l B®* T ^^--(A m x + B m r) - 
dx e 


e -U + 0- T ^ 

dx 


u ad — ©* T< f> ( x ) I = —A l BA®[ <f>i (x, r) + (p (x, r 


— Ida / 


(37) 


where (p (x. r ) is an approximation error which is to be kept small by a suitable choice of basis functions. 
Solving for u ad , we get 


M(ld — A 


dx 


A- { B®* t I (A m x + B m r ) + ©* T d> (x) — A _l fiA©7 <t > i (x, r) + (p (x, r) (38) 
dx £ ' 


From the assumption in Eq. (29), we can neglect the term ^©* T ^( A ) ) Then, one possible choice for the new 


basis function could be 


4-1 (*.«•) = f *(*) ^ 


(39) 


Alternatively, we can use the universal approximation theorem to approximate the uncertainty with a suitable 
choice of basis functions such as radial basis functions or sigmoidal basis functions 16 


dx 1 1 ~ 1 / T \ 

-y- = ~ (- A m x + B m r ) Bi©7d>i (x) + -B\8 (x, -) 

dX £ £ £ V £/ 


(40) 


where B\=eA BA and 

Since A m is Hurwitz and if ©, is bounded, then the Tikhonov’s theorem guarantees that the reduced solution with 
£ > 0 converge to the solution of the original system with £ = 0 as e — + 0. 


III. Optimal Control Modification Adaptive Law 

The tracking error equation in slow time is obtained as 

(x,t)-d(x,^j 


de dx„, dx 1 , 1 

—r = —7 ~r — —A m e H — B\ 

dx dx dx £ £ 


©[Oi 1 


We are interested in seeking an update law for © that minimizes the following cost function in slow time 

J= lim — [ (e — A) T Q(e — A)dx 
T /-“ 2e J o 


(41) 


(42) 


subject to Eq. (41) where A represents the unknown theoretical lower bound of the tracking error and Q = Q r >06 

R" x ". 

This optimal control problem can be formulated by the Pontryagin’s Maximum Principle. Defining a Hamiltonian 


H 


(e, ©}<I>i) = {e — A) T Q(e- A) + ^ p T (A„ 


e + Bi© 1 r< t , i -Bid 


2 £ ' ' ' £‘ 

where p (t) : [0,°°) — > R" is an adjoint variable, then the necessary condition is obtained as 

g = -v«,T = -i e( .- 


(43) 


( 44 ) 
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with the transversality condition p (xy) = 0 since e (0) is known. Treatng 0 <t> | as a control variable, then the 
optimality condition is obtained by 

(45) 


-V' 

A tp i £ 

The adaptive law which provides an optimal control solution can be formulated as a gradient update law as 


V// 07<t>, = ~P ' 


d® 


1, 


= -rvH & j = -r$i = ~r<t>i p J B { 

is an adaptive gain matrix. 


where T = T > 0 £ ' 


(46) 


The solution of p can be obtained using a “sweeping” method ' by letting p = Pe + S ® , <f>i , where P = P > 0 £ 
x " and S £ R' !X ". Substituting into the necessary condition yields 


dP 1 / ~ -r _\ C/kb ~r 

— — e H — P I A m e + B® j cfq — B\ S ) H — — 0j < t > i + S 
d T e V J d t 


c/5. 




c/t c/t e~' e 

This results in the following equations obtained by a method of separation of variables 


^Q(e-A)-^Al (Pe + S©;^) (47) 


f + I(M m+ A>) + I G -0 

( 

-fpPi ^©j T <I> + <5^) +5 


^ + - to+m)=o 


c/(©7<i>) i 

— p — QA = 0 

c/t e 


(48) 

(49) 

(50) 


For an infinite time-horizon problem when t y — > °°, then P (t) — > P (0) and 5 (t) — > 5 (0) for all t £ [0,°°). So, both 
P and 5 can be approximated by their constant solutions where 


PA m +A m P — — Q 
S = - A m J PB , 


(51) 

(52) 


Without loss of generality, a weighting constant v > 0 € R. is introduced to allow for adjustments of the modifica- 
tion term in the adaptive law. Then, v = 1 gives an optimal solution. Thus 


Then, the adjoint p is now expressed as 


S=-vA-'PB , 


p = Pe-vA m T PB i©7$i 


Substituting Eq. (54) into the gradient-based adaptive law yields the adaptive law in slow time 

H 5 = 4 r K T ™' - 


(53) 


(54) 


(55) 


Converting to regular time by multiplying e through Eq. (55) results in the optimal control modification adaptive 


law 


©i = -r 




PB\ — ©jP, PA m B\ 


(56) 


A. Stability Proof 

We now prove that the optimal control modification adaptive law (56) is stable and results in uniformly bounded 
tracking error. Toward that end, choose a Lyapunov candidate function 

V =e T Pe + trace (©[r^ 1 ©^ (57) 

Evaluating dV /d T in slow time yields 

^ = i e J [pA m +AlP^jeA^e J PB l (©[<f>i - s) - ^trace (©7<I> ie T PPi - vOIOiO^iPIpA^Pi) (58) 
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By the trace property 


trace 




= YX 


(59) 


where X,Y £ 
dV 1 


, then 
2 


z, 2 2 2 

= - — e T ge + -e J PB x ©7 $ - -e T PB! B - -efPB©^! + - v4>7 @iB7 PA~ '607O| 
dx £ e e e £ 


+ -v4>i r ©iB ] r PA m 1 B©[4> 1 (60) 

PA m 1 can be decomposed into a symmetric partM = 7 (PA“ 1 +A“ T P) = — iA“ T gA,„ 1 < 0 and an anti-symmetric 
part (V = 7 (PA” 1 — A m T P). Then, BA,,, 1 = M + N. By the property of a symmetric matrix, if the symmetric part of a 
matrix is negative definite, then the matrix is also negative definite. Since M < 0, therefore PA m 1 < 0. Thus 

dV 12 2 / i \ 2 

_ = __ e T ee _i e T rai5 + ^ v(J) T@ iB T ^--A“TgA“i +ivj Bi©^! + -V^eiB^BA'lB!©}^! (61) 

Using the property of an anti-symmetric matrix y : Ny = 0, dV /dx becomes 


dV_ 

dx 


^ = -ie T Qe-le T PB l 8-^v<t>J& l BjA- J QA- 1 B l &J® l + ^v<i>J&lBjPA- 1 B l &]<i> l 


2 

— 1 
£ 


(62) 


dV /dx is then bounded by 

^<-lA mi 4e)|| e || 2 + ^(B)|| e ||||B 1 ||^-lvA min (b}A“ t 2A“ 1 B 1 ) ||© 1 || 2 ||<I> 1 | 


2 

V (5 max 

£ 


(b, M-'Bi) ■ 


©s © m 


(63) 


where B) = sup, ||B|j, ©jj = sup, ||©j||, and A and <7 denote the eigenvalue and singular value, respectively. 
dV /dx can also be expressed as 

dV 1 

-p < —Kin (Q) INI Ikll -2 A max (P) IIBjII Bo] 
dx £ 


-lvA m; „(B[A- T eA- 1 B 1 )||© 1 ||||<I» 1 || 2 


©1 -2vcr mra (B[ J PA- 1 Bi)© 0 


(64) 


To show that the tracking error e and the weight variation © are bounded, we require dV /dx <0. Thus, it follows 


that 


2WP)|HH1B 1 ||Bq 
Klin ( Q ) 

2a max (BjPA m 'B l )&* 0 

i|| > ~7 

^ min [B\A m T QA m 7 l B\ 


k > 


l©ill> 


(65) 

( 66 ) 


Hence, there exists a compact set 'A where 


^ = J (^e &) ■ ||e|| < 2 ^' nax Ml l^ 1 II 1 1 1 1 < - G max{BjPA m l Bi)®* Q 


Klin ( Q ) 


Knin (b | A )n QA m B i 


(67) 


that contains the origin e = 0 and ©i = 0. 

Then dV /dx < 0 outside the compact set 7f. Thus, any trajectory e and ©i starting in 7? will remain in A"’ for all 
f. 18 Therefore, the compact set 7? is an invariant set. 10 Also , any trajectory e and ©i starting outside the compact set 
will approach the largest invariant set 7? as t — > °°. 18 It follows by the LaSalle’s Invariance Principle that e and ©i 
are uniformly bounded. Thus, the optimal control modification adaptive law is stable. 

From the Lyapunov stability analysis, it is noted that in the absence of persistent excitation, i.e., <f>i <!>, =0 for all t 
and if Bo = 0, then it can be shown by the Barbalat’s lemma that dV /d T is uniformly continuous and e — > 0 as t — » °°. 19 
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In contrast with the ct-modification adaptive law, even when the persistent excitation is removed, the tracking error 
does not tend to the origin. 20 

Since e and ©i are bounded, the unknown theoretical lower bound of the tracking error A at t = tt — > °° is also 
bounded by 


l|A|| < 


A maxjP) H#1 1 
A min ( Q ) 


j3 + So + 


evri 

Gtnin (A m ) _ 


( 68 ) 


where ©* $ </Jgl and 


d(@ T ‘ t>) 
dT 


< t] € M for all t. 


One unique feature of the optimal control modification is that as the adaptive gain increases and for V = 1, the 
system is guaranteed to be bounded. To show this, the optimal control modification adaptive law can be written as 


( 


©i = -rT>! [e' 


Then 


©! T>1 =-Bj (^Pe — vA„ l 1 PBi©^ r <J>i^ <J>i r<t>i 
Note that <t>7 r4>i € R, so for large <!>, r<J>i 


lim 


©7<I> 


l^i d t 


cf)^ r^>i — r<i>i 


= -B[ I Pe - vA m r PBi ©7 011=0 


.) 


the adaptive signal ©7 <t>i remains bounded and tends to 


Bi© 70 ! -^-p- l Aj n Pe 


The tracking error then becomes 


e -► ( A m + —P~ [ A, n P ) e — B\®\ Oj (x,t) - B\8 (x,t) 


which can be written as 


e-+-p- 


1 + V 

2v 


Q- 


1 — v 
2v 


e-B 


(©| 1 Oi (x,t) + S(x,t)^J 


(69) 

(70) 

(71) 

(72) 

(73) 

(74) 


where S =Aj / P — PA m . 

In a special case when 8 (x,t) = 0 and <I>i (x) = x, the system tends to a linear system as Oj TOi — >°°. Furthermore, 
if v = 1, then 

e(s) -+-H(s) B&fx(s) (75) 

where system transfer function matrix H (s) = (si + P~ l Q) * is strictly positive real (SPR) since H ( j (0 ) +H T (—/ft)) > 
0 as a result of P 1 Q > 0. For a SISO system, the Nyquist plot of a strictly stable transfer function for a SISO system 
is strictly in the right half plane with a phase shift less than or equal to j, 19 corresponding to a phase margin of at least 
f . Thus, one can deduce that the optimal control modification adaptive law is robustly stable with a large adaptive 
gain. 


B. Example 

Consider the following simple scalar system 


x = ax + bu + 6*x + v (t ) 


with actuator dynamics 

it = eA (u — u c ) 

where a < 0, A < 0, e > 0, |eA| <\a\, and v (?) is a disturbance signal. 
The reference model is 


Xm — U/nXm T bm^ 


(76) 

(77) 


( 78 ) 
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where a rn < 0. 

The actuator command is designed as 


Cl / Cl fit 

b\e A 


x + — & 1 <I> (jc, r) 

b £A 


(79) 


where <f> (x, r) = 


i T 


Note that if actuator dynamics are fast then the actuator command is 

bn 


-1 )x+—r-dx 
b V a / b 


The optimal control modification update law for slow actuator system is 


0 = -r (®epbi - vOT> T 0!^-^-j = -er (®ep^- - £V<M> T © 
where b\ = and p = — — , and for fast actuator system i 

ci 1 a m *' 


b 2 l 2 p 

\ ‘ 1 7 

\ a ci cipji 


(80) 


(81) 


e = -r 


xepb — Vx 2 db 2 — 


(82) 


If a and A are nominally in the same order of magnitude, then we note that for the slow actuator system, the 
effective adaptive gain is also reduced by £ for a similar performance as that for the fast actuator. 

For the numerical example, a = — 1, b = 1, 0* = 0.1, A = —1, £ : 0.1, a m = —5, b m = 1, r(t) = sinf, v(t) = 
0.05 sin lOf. The responses due to the standard MRAC adaptive law and optimal control modification adaptive law 
with the singular perturbation approach are plotted in Fig. 1. The response for the standard MRAC exhibits more 
initial transient than that for optimal control modification using the same adaptive gain. 



Fig. 1 - Responses due to MRAC and Optimal Control Modification 

Figure 2 is a plot of the control input and actuator command with the singular perturbation approach. As can be 
seen, the actuator command signal is quite large relative to the control input. This is due to the fact that the actuator 
dynamics are slow so a large actuator command does not translate into the same amount of control input for a finite 
time. The effectiveness of the optimal control modification is demonstrated by reducing the amplitude of oscillation 
in the control input significantly over that due to the standard MRAC. 
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40 



t, sec 



Fig. 2 - Control and Actuator Command due to MRAC and Optimal Control Modification 

Figure shows the responses due to the unmodified actuator command for fast actuator dynamics. As can be seen, 
the control input is insufficient to allow the plant to follow the reference model even with adaptive control. 




Fig. 3 - Responses due to Un-modified Actuator Command for Slow Dynamics 

IV. Flight Control Application 

Consider the following inner loop adaptive flight control architecture as shown in Fig. 1. The control architecture 
comprises: 1) a reference model that translates rate commands into desired acceleration commands, 2) a proportional- 
integral (PI) feedback control for rate stabilization and tracking, 3) a dynamic inversion controller that computes 
actuator commands using desired acceleration commands, 4) an adaptive controller with the conventional MRAC 
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law and with or without the optimal control modification adaptive law, and 5) a parameter estimator for actuator 
dynamics as relative to its nominal dynamics via the parameter e. The adaptive controller is designed to increase 
performance of the nominal dynamic inversion controller under adverse flight conditions such as upsets and damage. 
Under nominal fast actuator dynamics, both the dynamic inversion control and the adaptive law are computed without 
any differentiation between actuator commands and control inputs. As actuator dynamics degrade based on the value 
of e, both the dynamic inversion controller and the adaptive controller are modified accordingly to increase the actuator 
command signals. 



Fig. 4 - Adaptive Flight Control Architecture 
The linearized equations of motion are expressed as 

x = Anx+Ai 2 Z + Biu +fi (. x,z ) 


where A, 7 and B,, i = 1,2 are known, x = 

lT 


Z = A21X + A22Z + B2U+f2 (x, z) 
T 

P q r 


(83) 

(84) 


is a vector of roll, pitch, and yaw rates; 


z = 


A(j> A a A/3 AV Ah AO 


angle A/3, airspeed AV, altitude Ah, and pitch angle AO: u = 


is a vector of perturbation in the bank angle A 0, angle of attack A a, sideslip 

1 T 


A 8 a A 8 e A 8 r 


is a vector of additional aileron. 


elevator, and rudder deflections; and f) (x,z), i = 1,2 are parametric uncertainties which can be expressed as 


fi (• *,z) = ©* t <F(x,z) + § (x,z) 


(85) 


where 8 (x,z) is an approximation error which is assumed to be small by a suitable choice of basis functions, <f> = 
1 T 


Cl Ci C3 C4 


is a basis function for a sigma-pi neural network with Cj,i= 1 , ... ,4, as inputs consisting of 


control commands, sensor feedback, and bias terms; defined as follows: 

Ci = -p(h)V 2 x T ax T /3x T 


C2 = ^ p(li)V 2 y 1 </> 0 a fi a 2 /3 2 a/3 

C?, = -p{}i)V 2 u(x,z) T au(x,z) T /3m(x,z) t 


C, = -p{h)V^ 


px T qx T rx T 


( 86 ) 

(87) 

( 88 ) 
(89) 


where 0=0+ A0, a = a + Aa, /3 = /3 + A/3, V = V + AV, h = h + Ah, and 6 = 0+ A 6: and the overbar symbol 
denotes a trim state. 
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These inputs are designed to model the parametric uncertainty that exists in the damaged aircraft plant dynamics. 
For example, the aerodynamic force in the x-axis for an aircraft is given by 


pb 


- T + v q s ( Clq + C La a + C Lp P + C Lp ^ +C Lq ^ + Cl,. — + C Lga 8 a + C Lge 8 e + C Lgr 8, 


rb 

'TV 


a 


pb 


- ( Cd 0 + Co a a + Cd p ttt + Cd„ ttt. + Cd,. — + Co Sn 8 a + Cd g 8 e + Co Sr 8 , 


rb 


’2V 


<2V 


TV 


JSr Vr 


(90) 


Thus, Ci, C>, and C3 are designed to model the product terms of x, z, and u in the aerodynamic forces and moments 
equations; and C4 models the gyroscopic cross-coupling terms of x in the moment equations. 

The inner loop rate feedback control is designed to improve aircraft rate response characteristics such as the short 
period mode and the dutch roll mode. A second-order reference model is specified to provide desired handling qualities 
with good damping and natural frequency characteristics as follows: 


(s + 2t^p(OpS + (Op) <j) m — S p 8l at 
(s +2 i^qCOqS + 0) q ) 9 m = §q8[ on 
(^” -f- 2 ^ r CO r S T fit;) /)/;/ = §r8rud 


(91) 

(92) 

(93) 


where ® m , 0 m , and \f/, n are reference bank, pitch, and sideslip angles; co p , (O q , and ft), are the natural frequencies for 
desired handling qualities in the roll, pitch, and yaw axes; £ p , and t, r are the desired damping ratios; 8 /,,, , 5/ 0 „, and 
8 ru d are the lateral stick input, longitudinal stick input, and rudder pedal input; and g p , g q , and g r are input gains. 

Let p m = <f ) w , q m = 9 m , and r m = — j3,„ be the reference roll, pitch, and yaw rates. Then the reference model can be 
described as 


Xm — KpXm K( I Xmd T ~\~ Gr 


(94) 


where x m = 


Pm fit 

-I T 

8lat 8i on 8 rlu j 

Suppose, the elevator actuator is a slow actuator where 


, K p = diag (2 G) ; , , 2 t, q co q , 2 co , ) , K t = diag (a)|, co|, ft), 2 ), G = diag(g p ,g q ,g r ), and 


8 e = el e (8 e - 8 ec ) 


(95) 


where X e < 0 is the original elevator actuator rate, 8 ec is the elevator deflection command, and e > 0 is unknown 
but can be estimated by a suitable parameter estimation technique such as the recursive least-squares method if the 
control signal u {t )and its derivative it (t) is available by estimation, in which case u (t) is replaced by it (, t ), or by direct 
measurement. A recursive least-squares algorithm can be used to estimate e as follows: 


\ 


e = RAu I it 1 — Am t A t e 


R = —RAuAu J R 


(96) 

(97) 


where A u = u — u c . 

Assuming the pair (A 1 1 , B \ ) is controllable and z is stabilizable, the reduced-order equation for the pitch rate is 


0 = £2^1 1 A 1 1 x + EoB | l A\ 2 Z + E 2 U + E 2 B l l f\ 

are basis vectors. 


where £1 = 

1 0 0 

,£2 = 

0 1 0 

, and £3 = 

0 0 1 

Note that E 211 = 8,. and £2M e = 

= S ec , then differentiating the reduced 


£2^1 1 A 1 ix + £ 2 £ 1 l A n z = eX e (E 2 B l l A n x + E 2 B l l A l2 z + E 2 B l 1 /i+£ 2 M f ) 
Setting x = — ^ K p + x + Gr, then the actuator commands can be computed as follows 


(98) 


( 99 ) 


U c = 

£i£r ] 

A;E 2 Bi 1 A n 


— f K p H — -jx + Gr 

+ 

0 

TT c E 2 B \ 1/1 12 

z-B 1 ^Ahx+Ai 2 Z + ©i 


£3V 

3x3 

\ S / 


0 

3x3 
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where u c = y A5 ac A S ec AS rc 

The closed-loop angular rate dynamics are then given by 


x = — I K„ H — - )x + Gr — 


E i 
£ 3 


(©[«>- 5i) 


( 101 ) 


3x3 


Let e = 


So x)dz x m x 


be the tracking error, then the tracking error equation is given by 


e = A m e + B — 8\ 


where 


Ani 


0 / 
-K, —K„ 


B = 


0 

E 


E = 


Let Q = 21, then the solution of Eq. (51) yields 


Ei 

eX e E 2 A j"/ 
£3 


3x3 


P = 


Kr'K p + K-' (Ki + I) Kr 1 

K7 


K p l (l + Kr') 


>0 


A m l is computed to be 


AZ = 


- Kr l K p -Kr 1 


Evaluating the term B r PA m l B yields 


B T PA7„ l B = - E T Kr 2 E < 0 


( 102 ) 

(103) 

(104) 

(105) 


Applying the adaptive optimal control modification (56), the weight update law is then given by 


©1 = -r <i>(e T PB + v<i> T ® l E T Kr 2 E 


(106) 


(107) 


(108) 


(109) 


A. Simulation Results 

A simulation study was conducted using a generic transport model (GTM) which represents a notational twin-engine 
transport aircraft as shown in Fig. 5. 21 An aerodynamic model of the damaged aircraft is created using a vortex lattice 
method to estimate aerodynamic coefficients, and stability and control derivatives. For the simulation, a damage 
configuration is modeled corresponding to a 28% loss of the left wing. The damage causes an estimated C.G. shift 
mostly along the pitch axis with Ay = 0.0388c and an estimated mass loss of 1.2%. The principal moment of inertia 
about the roll axis is reduced by 12%, while changes in the inertia values in the other two axes are not as significant. 
Since the damaged aircraft is asymmetric, the inertia tensor has all six non-zero elements. This means that all the 
three roll, pitch, and yaw axes are coupled together throughout the flight regime. In addition, the elevator actuator is 
simulated as an impaired flight control actuator with a much smaller bandwidth than that for a healthy actuator. 
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Fig. 5 - Left Wing Damaged Generic Transport Model 

A level flight condition of Mach 0.6 at 15,000 ft is selected. Upon damage, the aircraft is re-trimmed with T = 
13,951 lb, a = 5.86°, 0 = —3.16°, 8 a = 27.32°, 8 e = —0.53°, 8, = —1.26°. The remaining right aileron is the 
only roll control effector available. In practice, some aircraft can control a roll motion with spoilers, which are not 
modeled in this study. The reference model is specified by co p = 2.0 rad/sec, co q = 1.5 rad/sec, ft),. = 1.0 rad/sec, and 
Cp = C<? = Cr = 1 / \/2. The simulations also include a random signal to represent sensor noise. 

The actuator dynamics are modeled with X a = K = X r = 50 /sec with position limits of ±35° for the aileron and 
elevator and ±10° for the rudder. The eparameter for the elevator is set to 0.01. 

The pilot pitch rate command is simulated with a series of ramp input longitudinal stick command doublets, 
corresponding to the reference pitch angle ±3.81° from trim. At t = 10 sec, a wing damage and elevator actuator 
degradation occur. The tracking performance of the baseline flight controller, which is a proportional-integral feedback 
type with no adaptation, is compared against the optimal control modification adaptive law with and without the 
singular perturbation approach for slow actuators. Both the adaptive laws are implemented as an augmentation to the 
baseline controller. An adaptive gain of F = 60 and a weighting factor v = 0.2 are selected. 

The aircraft angular rate responses are shown in Figs. 6 to 9. Figure 6 illustrates the pitch rate responses. With 
no adaptation, the baseline controller cannot follow the reference pitch rate very well. The pitch rate response clearly 
lags the reference model significantly due to the degraded elevator actuator. The optimal control modification with and 
without the singular perturbation approach significantly improve the tracking and reduces the lag between the response 
and the reference model. The pitch rate is worse without than with the singular perturbation approach, as large initial 
transients occur at failure. So, the singular perturbation approach demonstrates an improved response due to the slow 
elevator actuator. 

Since the damage occurs to one of the wings, the roll axis is most affected. With no adaptation, there is a significant 
roll rate as high as 20°/sec as shown in Fig. 7. There is a steady-state oscillation of ±15°/sec in the roll rate. Both the 
optimal control modification adaptive laws reduce the oscillation to within ±2°/sec. However, there is a large initial 
transient of about 20°/sec due to the sudden wing loss. 

Figure 8 is the yaw rate response of the damaged aircraft. The optimal control modification with the singular 
perturbation approach slightly improves the yaw rate response than that without. Both adaptive laws result in a much 
lower yaw rate than the baseline controller. 

Figure 9 is the plot of the tracking error Jz ?2 norms for the three axes in combination. Not surprisingly, without the 
adaptation, the baseline controller suffers a large tracking error. 
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Fig. 6 - Pitch Rate 



Fig. 7 - Roll Rate 
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Fig. 8 - Yaw Rate 





Fig. 9 - Tracking Error Norm 

The control surface deflections are plotted in Figs. 10 to 12. The aileron actuator is a fast acting actuator, so in 
all cases, the aileron deflection tracks the command as shown in Fig. 10. With both the optimal control modification 
adaptive laws, there is a large aileron command at the instance of failure that causes the aileron to saturate. The spike 
in the aileron command is larger without than with the singular perturbation approach. 


16 of 21 


American Institute of Aeronautics and Astronautics 


40 




Fig. 10 - Aileron Deflection 

Figure 1 1 is a plot of the elevator deflection. Since the aileron actuator is degraded, it is clear that the elevator 
command cannot be tracked well by the elevator actuator. The optimal control modification adaptive law without the 
singular perturbation approach produces a large initial spike in the elevator command of about —25° while the initial 
transient in the command is reasonable with the singular perturbation approach. 



Fig. 1 1 - Elevator Deflection 

The rudder deflection is shown in Fig. 12. With no adaptation, the rudder deflection oscillates from about 0° 
to —7°. Both the optimal control modification adaptive laws produce very similar rudder deflections, which reduce 
steadily from the initial transient to its new trim value of —1°. 
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Fig. 12 - Rudder Deflection 

The attitude responses of the damaged aircraft are shown in Figs. 13 to 16. When there is no adaptation, the pitch 
attitude could not be followed accurately as seen in Fig. 13. With the adaptation on, the tracking is much improved and 
the optimal control modification adaptive law performs better with than without the singular perturbation approach. 





Fig. 13 - Pitch Angle 

Figure 14 is the plot of the bank angle. Without the adaptation, the damaged aircraft exhibits a rather severe roll 
behavior with the bank angle ranging from —40° to 20°. With the adaptation on, the roll angle is drastically reduced 
from an initial transient of about —30° without and —25° with the singular perturbation approach to almost zero. 
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Fig. 14 - Bank Angle 

Figure 15 is a plot of the angle of attack. The baseline controller results in a significantly large initial angle of attack 
when the damage occurs. The maximum angle of attack is about 12°, which could be close to stall. The maximum 
angle of attack is reduced to 6° when the adaptation is on. 




Fig. 15 - Angle of Attack 

Figure 16 shows a plot of the sideslip angle. The baseline controller produces a steady oscillation in the sideslip 
angle between ±2°. With both the optimal control modification adaptive laws, the sideslip angle is reduced to near 
zero. 
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Fig. 16 - Sideslip Angle 

V. Conclusions 

This paper presents a singular perturbation approach in connection with an optimal control modification adaptive 
law for a time-scale separated system with slow actuator dynamics. The singular perturbation approach transforms 
the system into a reduced-order system in a slow time coordinate. The actuator command is obtained by the model 
matching condition in the slow time coordinate. The resulting actuator signal in effect is increased by the ratio of 
the norm of the plant’s transition matrix to the norm of the slow actuator’s transition matrix. The optimal control 
modification adaptive is derived and analyzed for stability using the Lyapunov method. Under fast adaptation when the 
adaptive gain is large, the analysis shows that the tracking error remains bounded and stable. The singular perturbation 
approach with the optimal control modification adaptive law is extended to a flight control application. Simulations of 
a flight control performance for a damaged aircraft with an impaired elevator actuator demonstrates the effectiveness 
of the method. 
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